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SllviMARY 

A method I 3 presented for ehe calcuj,at;on of elastic stresses 
In syinn/e tr i cal disks typical of those of a high -temperature gas 
turlLne, The method is essentially a f Inite-difference solution of 
the equllihrium and compat ibility equations for elastic stresses in 
a s;>'mnetri cal disk. Account can be taken of point -to -point varia- 
tions In disk thickness, in temperature, in elastic modulus, in 
coefficient of thermal expansion, in material density, and in 
Poisson’s ratio. No numerical integration or trial -and -error pro- 
cedures are involved and the computations can be performed in rapid 
and routine fashion by nontechnical computers vith little engineering 
super\isLon, Checks on problems for which exact mathematical solu- 
tions are known indicate that the method y'elds results of high 
accuracy . 


Illustrative examples are presented to show the mariner of 
treating solid disks, disks with central holes, and disks con- 
structed either of a single material or of two or more welded 
materials. The eifect of shrink fitting is taken into account by 
a very simple device. 


INTRODUCTION 

One of the problems in the design of gas turbines is the deter- 
mination of the stresses in the turbine disk under operating con- 
ditions. Calculation of the elastic-stress distribution is a first 
step in the determination of the true stress distribution. This 
stress distribution is based on the assumption of linearity of stress 
with strain and differs from the true stress distributions, which 
may contain stresses beyond the proportional elastic limit of the 
material . 

The eqiiations lor the ela,stic-atross distribution in symmetrical 
disks are well known. Their solution may, however, offer considerable 
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b rim of disk or base of blades 

Example of the use of double subscript; 

redial stress Oj. at station (n-l) 

The following supplementary s;ymbol 3 denote combinations of the 
foregoing symbols arising in the analysis: 
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AssuzaptionB 

Tbo assumptions are made tJiat stress Is prcporticnal to si-rain 
and that the disk material is complete l;v' elastic at the stress die- 
trllution induced ty the centrifugal and thermal effects. All 
variahles of material properties and operating conditions are assmed 
to be symmetrical about the axis of rotation. Axial stresses are 
neglected and at any radius, the radial and tangential stresses are 
assumed to be uniform across the thickness of uhe disk. Temperatures 
are taken in the central piano perpendicular to the axis of the disk. 


Outline of Method 

In a thin rotating disk of variable thickness, the state of stress 
at any radius can be completely defined by the two principal stresses, 
the radial and tangential stresses O-^ and Ot • equations are 

therefore necessary to determine the two unknown stresses. The first 
of those equations can be obtained from the conditions of equilibrium 
of an element of the disk; the second, fi‘om the compatibility con- 
ditions, vrhlch are mathematical statements of the interrelation between 
the radial and tangential strains in a symmetrical disk. 

The equilibrium and compat ibility equations result in differen- 
tial form defining relations between the stresses at radius r and 
those at a radius inf initeslmally removed from r. Except for some 
special cases, the solutions of these equations are difficult to obtain. 
In order to facilitate solution, the differential equations are rewrit- 
ten in finite difference form relating the stresses at radius r with 
those at a radius finitely removed from r. By means of the finite- 
difference equations, the stresses at an arbitrary finite number of 
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Btatlone along the disk radius are expressed in terms of the stresses 
at a single reference station near the center of the disk. For a 
disk with a central hole the reference station is chosen at the inside 
radius, where the radial stress is zero; hence, the stresses at all 
stations in the disk are expressed in terms of the single unknown, 
the tangential stress at this station. For a solid disk, the refer- 
ence station is chosen at a point near the center of the disk (at a 
radius of about 5 percent of the disk radius). In this region the 
radial and tangential stresses can he assumed to he approximate l^^ 
equal; again, therefore, the stresses at all stations are expressed 
in terms of a single unknown. 'The unknown can then he determined hy 
the boundary conditions at the rim of the disk where the radial stress 
is equal to the centrifugal bucket loading. When the radial stress 
at the rim, expressed in terms of the tangential stress at the refer- 
ence station, is equated to the bucket loading, the tangential stress 
at the reference station is evaluated. After the tangential stress 
at the reference station has been determined, all the other stresses, 
expressed in terms of this stress, can be evaluated. 


Differential Equations 

The equilibrium equation as given in reference 2, (p. 374) using 
the notation of this paper, Is 

~ (rh Qr) - h at + pCi)^ = 0 (l) 

The compatibility relation is obtained by elimination of u 
(the radial displacement of any point, on the disk as the disk passes 
from the unstressed to the stressed condition) from the stress-strain 
displacement equations 
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Equation (3) is subtracted from equation (2) to eliminate u, 

du _ u ( 1 + M-) (O'r ~ ) 

dr r 
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( 4 ) 
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But 

(S) 


(7) 

Equations (l) and (7), together with a knowledge of the boundary 
conditions, are sur'ficlent to solve for the two unknowns O^. and 
*^t • ^Because p, q, a, AT, and h are, in general, functions 

01 the radius r, the equations cannot readily he solved in their 
difxerenbiai. tom; a fiiiits-difference solution was therefore derTved . 


rdu — udr _ d / u 
r 2 dr ~ ^ \r / 


I’herefore, by equations (3), (5), and ( 6 ), 

d /®t\ d /h CT}-\ 5 (i 

dr \ E / dr dr 


Finite -F.l.i'f erence Equations 

Tile trsansla-ion oi dl '.TerenLial equations into f inite-difrerenco 
form to faoiii t.ate nciutiori la common in engineering practice . The 
metnod liaa. in fac'c^ been apolieci in limited fashion to the solution 
of the steam- lu.rblne disk: problem (reference 2, pp. 598-4-00). This 
application neglects^ however^ the poirit-to^-pcint variation :ln 
phjsica.*. prepertioG, and therefore no application to the gas-turbine 
dislc in wiiich there 1 vS appreciable variation in properties from hub 
to rim is ?nade. In addition^ the solution of the equations involves 
an interpolat;Lon procedure the elimination of which could reduce 
con3idex‘'aD.^.,y the ams;unt of calculation necessary for a solution and 
increase the .'iCCL.racy of the final results, 

A number of discrete point stations are chosen along the disk 
radius as shoim In figure 1(a). If ft is assumed that the stress 
distribution In the disk has alree.dy been determined^ all quantities 
equations ( 1 ) and (?) are therfore known at each of the 
point stations and ohe values o.f corresponding quantities at the 
point A midvay bo ween the nth and (n-l}st point stations can 
ohen be appro xirnately dcteimined * For example^ in the plot of rh 

against r (fig. 1(b)), the radius at point A is expressed as 

n = I + On) 
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the value cf rh is 

(rh ^n-1 ^r,n-i -n ^r,n^ 

and the slope cf the curve at point A, which is approximately equal 
to the slope of the chord ioining points n and n-lj is defined as 
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In a similar way the values of each of tne other variables entering 
into the equations can he evaluated for point A. If the evaluations 
are correct^ the quantities at A must satisfy equations (l) and (7), 
These equations therefore hecomo, in finite -difference fona 
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Solution of Finite -Difference Equations 

Equations (lO) and (ll) represent two equations from which 
^r n ^t n expressed in terms oi' ^^n-1 ^t,n-l* 

If the linear nature of the equations and the posslhility of 
successive application of the equations to proceed from one 
station to the next are considered, the stresses at any station 
can ultimately he expressed in linear terms of the stresses at any 
other station. It will be convenient to express the stresses at 
all stations in terms of the stresses at the station a. At this 
station, the unknown value is the tangential stress hence, 

the stresses at station n are expressed in the linear^ terms ^ 

^r,n ~ *^r,n ^t,a ^r,n 

^t,n ' ^t^n ^t,a ®t^n • (12a) 

and those at station n-1 in the form 

^r,n-l “ ^r,n-l ^t,a ^ ®r,n-l 
^t,n-l = H,n-1 ^t,a + ®t,n-l 


where the coefficients A^^^n; ®r,n; ^t,n 

to be determined. 


The substitutvon of equations (l2a) and (ifb) into equa- 
tions (lO) and (ll) and the separation of the terms with and without 
^t,a equations 


(Cn Ap^j^ D^ ^t,n ^r,n-l ^n ^t,n-l) ^t,a 

+ (^n ®r,n " ^t,n ^ ^n ®r,n-l “ ^n ®t,n-l ^ = 0 (13) 


and 


(o’n ^r,n ^’n ^r,n-l ^’n ^t,a 

+ (C'n ^r,n ^ ^'n ^t,n “ ®r,n-l ^ xi ®t,n-l ~ ^*n^ “ ^ 


(14) 


The stress ^ is really arbitrary as far as equations (15) 
and (14) are concerned because it depends upon the boundary conditions 
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and not on the equations of elasticity from vhich equations (13) 
and ( 14 ) were derived; that hy a suitable choice of the factors 
that determine boundary conditions, such as bucket loading and 
shrink fit, ^ can be set at any desired value without Inval- 
idate r:ig in ai\y way the equations of elasticity (l) and (7), or their 
ultimate finite-difference forms in equations (13) and (14). If an 
equation in the form cx + d = 0 is to be true for all values of 
X, the coefficients c and d both must be zero. Because equa- 
tions (13) and (14) are to be true Independent of the value of 
the coefficients of ^ must be zero, and the two equations 
reduce to ^ 


-^r,n “ -^t,n ^r^n-1 ‘^,n-l ” 

^’n -^t,n ~ ^’n ^,n-l ^t,n-l “ ^ 
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C’n Sr,n - “’n »i,n Bj-.n-l + O'n Bt.n-l " H’n « 0 


from which A^^n; ^r n; Bt^u can be determined in the 

f orm 
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(16) 


If the c ’efficients A^ ®r n? B-|-^^ are known 

for station n-1 they can be determined by means of equation (16) 
for station n. 


The coefficients at the first station (r = a) can be determined 
by inspection for both the solid disk and the disk with a central hole. 
Inspection of equation (12a) shows that for a solid disk in which both 
the tangential and radial stresses at the first station are equal to 
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For disk with a central hole in which the radial stress at the 

first station is zero and the tangential stress is 0+ _ 

a 


•^r,a ~ ®r,a 


®t,a = 0 


■^t,a “ ^ 

From these known coefficients at the first station, the coefficients 
at all other stations can therefore ho determined bj successive 
applications of equation (16). Once all the coefficients have been 
determined, the unknown can be determined. The radial stress 

at the rim is the centrifugal loading of the buckets 

*^r,b “ ^,b ^t,a ^ ®r,b 


or 
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t,a 


^r,b ^r,b 

Ar,b 


(17) 


where ®r,b coefficients for radial stress at 

the rim. The radial and tangential stresses at all stations can be 
obtained from equation (l2a) after ^ and all the coefficients 
have been determined ’ 


ILLUSTRATIVE APPLICATIONS 

Case I - Elastic-Stress Distribution in Solid Disk 

The profile of a disk that is to be analyzed for stress dis- 
tribution at a speed of 11,500 rpm and the temperature distribution 
are shown in figures 2(a) and 2(b), respectively. The first step 
in the analysis is to choose an arbitrary number of stations along 
the disk radius. The first station is chosen at a radius of about 
5 percent of the rim radius, the last at the rim. The stations need 
not be equidistant, in fact, it is advisable to choose the stations 
closely together where there is sharp change in disk contour, in 
temperature graclent, or in variation of physical properties. In 
this case 18 stations were chosen, spaced relatively close together 
near the rim where the gradients in temperature and physical prop- 
erties were high and near the center for subsequent use of the same 
example to illustrate the effect of a central hole. When only a 
Solid disk is considered, no concentration of points near the center 
la necessary. Tiie various steps of the calculation are tabulated in 
table I. 
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The disk radius at each station is listed in column 1 of table I. 
The thickness of the disk h at each station is listed in crlumn 2. 

A sharp discontinuity in thickness, such as an abrupt flange, should 
be faired in the disk contour ana the faired disk used in determining 
thickness . 

Ordinarily the density of the material is constant throughout 
the disk, even over the wide range of temperatures. If a faired disk 
has been used, however, the density of the material in the faired 
region should be adjusted to produce the total mass that actually 
exists in the region of each station. Although a flange does not 
reduce the stress at its own region by increasing the area, its mass 
must be included as it produces centrif-ugal stresses throughout the 
disk. The corrected density at each station multiplied by the square 
of the rotational speed is listed in column 3. In this case no 
fairing was necessary; hence, all values of denniuy are equal. 

Poisson’s ratio, listed in column 4, has only an insignificant 
effect on the stress distribution and, because no accurate data are 
available, a constant value of 0.4 may be used. If accurate data on 
the variation with temperature on Poisson’s ratio are available, use 
of the exact variable values presents no greater difficulty than use 
of a constarit value. The values of [x used in this example are 
shown in figure 2(c), and were for convenience obtained by the assump- 
tion of a linear variation in \x with temperature. 

The modulus cf elasticity at each station is listed in column 5. 
Variations in this property have a significant effect on the final 
stress values and accurate data should be used if available-. For 
this example E was arbitrarily assumed to depend linearly upon the 
temperauure, and the variation along the radius is shown in fig- 
ure 2(d). In practical computations, bhe true values of elastic 
modulus associated with the particular temperature at each station 
may be used . 

The coefficients of thermal expansion are tabulated in column 6. 
These coefficients mu.st be the average values applicable to the range 
between the temperatures actually existing and those at which there 
is no thermal stress. For a homogeneous disk in which there is no 
shrink fitting of one part to another, the condition of zero thermal 
stress is at room temperature. Engineering tables usuall^f list the 
average temperature coefficient of expansion between room temperature 
and values of high temperature; the listed values may therefore be 
used directly. 

The difference between the actual temperature and the temperature 
at which there is no thermal stress is listed in column 7. In this 
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case I'.he stress-free condition is at a room temperature of 70^ F. 
Column 7 is theref<ore obtained hj eulDtracting 70^ F from each of the 
values in f igure 2(b). This column is of great significance in the 
case involving shrink; fits. 

The quantities to are computed for each station as 

indicated in columr.s n to 34 of table I. Values in each of these 
col.umns can be obtained in one set of oxjeratlons on a standard 
computing machine. The method of obtaining the data from the 
suitable previous columns is indicated at the heading of each column. 

Values In columns 33 and 34 must be computed simultaneously. 

The first value for each of those columns is unity. Subsequent 
values make use of the previously obtained values in the same columns. 
Thus, to determine the value for column 33 at station 2, column 27 
at station 2 is multiplied by column 33 at station 1, and the product 
is added to the product of column 20 at station 2 by column 34 at 
station 1. For exampl.e ; 

0.81902 X 1.0 + 0.10098 X 1.0 = 1.00000 


Columns 35 and 36 are likevlse computed simultaneously . The 
first value in each of these columns is zero and each subsequent value 
is obtained from the previous values in accordance with the symbolic 
notation given at the head of each column. Thus, to obtain the 
column 35 at station 2, column 27 at station 2 is multiplied by 

column 35 at station 1, column 28 at station 2 is multiplied by 

column 36 at station 1, and the two products are then added to 

column 31 at station 2 


0,81902 X 0 + 0.18098 X 0 - 67.192 = -67.192 


Column 37 is uniform for all stations and is obtained from the 
expression 

^r,b - 

(3VFX 


where Op bucket loading at the rim. The bucket loading is 

obtained by dividing the total centrifugal force at the root of the 
buckets by the total rim peripheral area. In this problem is 

6500 pounds per square inch, column 37 is therefore, 


8530 - (-72,896) 
i: 93765 


42,003 pounds per square inch 
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Columns 38 and 39 f ive tie radial and tangential stresses at eacli. 
of the stations. As indicated i.r'. table I they are obtained by routine 
iriultipll cat ions and additions of columiivS 33 to , 

The radial and tangential stresses from columns 38 and 38 are 
plotted in figure 3. The stresses at the center of the dislc are taken 
enual to those at station a, which io one-ljalf inch removed from the 
center. 


Because the method presenteO. is the only one knovm to the author 
that takes into account point-to-point v 'riation in Poisson’s r-itio, 
the error invo'^.ved in the assuiaption of a constan-t oalue of this 
quantity as compared with the rigorous troatoent c.:f its point-to-point 
variations is \al.uab1.e to determine. The brokea-l,:‘ne ciaryes in fig- 
ure 3 show calculations for oonsiant values of y 0,3 and \x - 0.5 
compared with the soiJ.d cuiuos^ \dilch ahov; the atre'^sea i“or a con- 
tj.nuously variable value of \i with tompei'oturo^ as slun^i in figure 2 
and tobulatel in coinrin 4 of taolo The near- coincidence of those 
curves indicates that toe assumqt lory of a co ist'^nt value of p vriLhiii 
the range of actual \alues results in accurate final vaj.uos of radial 
and tangential stresses. 

The effect that di'^ference i-u the nimler stations has on tl'.e 
accuracy of the results is snowr in figure 4, Little accuracy ±3 
g'^ined by tie use of additional points; as few six pc'ints in this 
particular case can yield accurate results at a gi'Cat neivlng in corn- 
put in^ time. 


Case II - Blast ic-Gtress Distribution in 
Disk with Centre 1 Hole 

A disk with a central hole is studied in a manner similar to tlie 
solid disk except that the first station is t^iion ai bhe inside loui-d- 
ary instead of at an arbitr-'ory small d i.stanoo as In the solid disk. 

The jholCG of stitioo.:: near bue central liole is, ho'-ever, critical for 
tins case. Stations should be taken also at dj ybarijes of 1, 2, 5, and 
5 porcont of the rim diemoter from the inoido boui-d^iry of the dink. 

In order to illustrate the procedure, the dick of 'L.ruro 1 is again 
used but with a conti/al hole 1 i.nch ^ri (l.iaineber, chooon so that 
station a will be ccnveuiently located in ui:3 name o^^ace as sta- 
tion a for Case I, 

Columns 1 to 32 for the dis:: of figure 2(a) with the central iiolo 
are identical to the correspondi:ig coliimns of table 1 for the solid 



NACA TW No. 1279 


15 


disk. In column 33 the entry for station a is 0 instead of 1 as 
for the solid disk; otherwise, the procedure for calculating col- 
umns 33 to 39 is the same as tliat of table I. Table II gives the 
modified columns 33 to 39 that result from changing the single 
first entry in column 33 and figure 5 shows a plot of the resulting 
radial and tangential stresses; the curves marked ”18 stations” are 
the stress values for this computation. 

The results of supplementary calculations using different 
numbers of stations (fig. 5) indicate that considerable error can 
result in the determination of the peak stress at the inner boundary 
if an insufficient number of stations are chosen near this boundary, 

A more , judicious choice of stations for the 6- and 10-station systems 
could produce more accurate results than those shown; in the absence 
of experience in choice of locations, however, it is better to choose 
a large number of stations and insure accuracy. 

A practical procedure used to reduce the amount of calculation 
necessary to obtain the critical end stresses is to calculate the 
stress distribution on the basis of a solid disk using a few stations 
and then to modify the stresses in the immediate vicinity of the 
central hole by the stress concentration factor characteristically 
introduced by the hole. A comparison of figures 4 and 5 indicates 
that, with the exception of the region Immediately adjacent to the 
central hole, the stresses are similar for the cases of the solid 
disk and of the disk with the central hole; this stress distribution 
depends very little on the number of stations chosen. By reference 3 
(fig. 145), for example, the characteristic stress concentration for 
a disk with a central hole of which the diameter is one -twentieth of 
the outside diameter is about 2.0. From calculations based on dif- 
ferent numbers of stations the calculated average stress at the center 
of the solid disk (fig. 4) is 43,000 pounds per square inch. The 
tangential stress at the inside boundary for the disk with the central 
hole should therefore be 2 x 43,000 = 86,000 pounds per square inch. 
The radial stress at a free boundary is, of course, zero. A curve 
faired between these boundary values and the general curves of fig- 
ure 4 would coincide very closely with the 27 -station result of fig- 
ure 5. 


Case III - Elastic-Stress Distribution in 
Composite Welded Disks 

For some applications turbine wheels must be fabricated by 
welding parts composed of several mater ial.s. The method of analysis 



16 


NACA TN No. 1279 


presented is applicable to studies of composite welded disks in which 
various alternatives of boundary location and shrink interference can 
be investigated with few changes in the tabulated computations. The 
procedure is Illustrated for a typical application in which the bound- 
ary location is constant. 

The disk of i Igure 2(a) Is assumed to be made in two parts with 
the boundary at the 6 -inch station. Figure 6 shows the two portions 
of the disk Just before welding. The heat-resisting outer portion 
is heated to 670° F while the inner portion is maintained at 70° F. 

(in practice both portions may be heated while maintaining a desired 
temperature differential.) At this temperature condition an exact fit 
exists between the mating tips of the two parts. The wedge is then 
filled with weld metal. 

The assumption is made that this temperature differential between 
the two portions of the disk is maintained throughout the welding 
process in making the calculations. Any cooling of the outer region 
prior to the placement of the weld metal would produce a crushing of 
the mating tips and reduce the effective amount of shrink. Localized 
effects of the weld metal in producing residual stresses are neglected. 
The calculations are made as if the high-temperature alloy, having 
full width at the mating face, is shrunk at 670° F onto the full-width 
steel central portion. Table III shows the essential tabulations for 
this case. 

In order to insure accuracy a few more stations than were used 
in tables I and II have been chosen in the vicinity of the boundary. 
The densities of the two materials are somewhat different. The 
quantities for p and E are the values at the temperatures of 
figure 2(b). The quantity AT at each station is the difference 
between the existing temperature and the temperature at which there 
is zero thermal stress. The temperatures of zero thermal stress occur 
Just before the shrink fit when the outer portion is at 670° F and 
the inner portion is at 70° F. Therefore, for the outer portion 
670° F is subtracted at each station from the temperatures of fig- 
ure 2(b) and for the inner portion 70° F is subtracted. The value of 
a at each station must be the average a between the stress -free 
temperature and the operating temperature. At the rim, for example, 
an average coefficient of expansion between 670° F and 1270 F must 
be used. The average temperature coefficient applicable to the 

range between any two temperatures Tp and Tg can be found by 


a*2 T'2 -a'l T'p 



°l-2 = 



NACA TN No. 1279 


17 


where average coefficients of thermal expan- 
sion between between room temperature and the temperatures and 

T2^ respectively^ and T']_ and T’^ Q-i'e the temperature differences 
between T]_ and Tg and room temperature . 

The procedure of calculating table III from column 8 on is 
similar to that of table I. The final calculated values of stress 
are shown in figure 6. 

Comparison of figures 5 and 6 shows tiiat a shrink fit^ unless 
excessive, can have beneficial effects. The shrink fit reduces the 
tangential tensile stresses that exist near the central hole under 
operating conditions and also the tangential compressive stress at 
the rim. Compressive stresses at the rim can be detrimental. If the 
elastic compressive stresses exceed the yield strength of the material, 
plastic flow takes place and a residual tangential tensile stress 
exists after operation. Because the region of the rim is a stress- 
concentrated area as a result of the blade attachments, even rela- 
tively small residual tensile stress may cause cracks. The shrink 
fit removes the high tensile stress at the center and the high com- 
pressive stress at the rim but Introduces a high tensile stress at 
the boundary of the two fitted regions. The boundary is at a lower 
runniiig temperature and has no stress -concentrating effects of the 
blade attachments. 

Ihe optimum amount of shrink, however, is I'airly critical. 

Probably the shrink of the illustrative example is excessive. There 
is no noed to reduce the stresses at the hole and at the rim as much 
as shown in figure 6 at tae expense of such a high stress at the 
boundary. An additional calculation can readily be made by using a 
smaller temperature difference of shrinking than assumed in this 
calculation. Only columns 6, 7, 24, 31, 32, and 35 to 39 are affected 
by any change, and the redistribution of stress can be calculated vei^ 
rapidly. Thus, a more suitable shrink fit can readily be found. 


Case IV - Check on Adequacy of Method 

Checks on the adequacy of the method were obtained by comparing 
the results of finite-difference calculations to theoretically correct 
results in several cases where the latter could be obtained. In one 
case a parallel sided disk was studied. The conditions of operation 
are shown in figure 7; the circles and squares show the radial and 
tangential stresses as determined by the finite -difference method, 
1‘es^ectiVely ; and the solid lines show the theoretically correct 
stresses obtained by rigorous solution of equations (l) and (?) for 
this case which E, a, and p are constant. This correlation 
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is seen to "be very good. The maximrin dev5ation occurs at the boundary 
of the central hole where the difference between the tangential stress 
as computed by the finite-difference method and the theoretically cor- 
rect value is about 2 percent. The average deviation between the 
theoretical and finite difference stresses throughout the disk is less 
than ione-half percent. Checks on solid disk produced closer agree- 
ment, even when a small number of stations were used. A check on a 
cisk of uniform strength produced results differing from the exact 
solution in the order of ±one-fourth percent of the theoretical 
stresses throughout the entire disk. 


CONCLUSIONS 

The finite-difference method of calculating stresses in rotating 
disks has been applied extensively to various types of turbine disks 
under different conditions of constant temperature or with a temper- 
attire gradient. The procedure was found to be convenient and rapid. 
Wliere checks were available, the results showed a high degree of 
accuracy. 


Flight Propulsion Eesearch Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, February 27, 1947. 
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TABLE I - CALCULATION OP STRESSES 

[Engine speed, 11,500 rpm; operating 
at center; average radial stress 
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table III - CALCtlLATIOH OF 

[Engine apeedi 11|600 rpm; operating 
at center; shrinking condition, 6 ’tCo F for 
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Figure !• 


lb) Value of typical function mldv^ay between 
point stations n and n-| , 

— Sketches used to derive finite— difference equations for 
stresses In symmetrical rotating disk. 
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Figure 3. — Stress distribution in solid disk of figure 2 

computed with constant and with variable values of Poisson's 
rat i o. 
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Figure 4. - Effect on calculated stress distribution in disk 
of figure 2 of various numbers of stations. 
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Figure 5. * Effect on calculated stress distribution In disk of figure 2 with 

central hole of various numbers of stations, (Location of stations shown 
on fig, 4,) 
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Figure 6. - Stresses at running conditions of speed and temperature in com- 
posite welded disk. Temperature of operating rim section during welding, 
670® F; temperature of central section, 70® F. 



Stress, Ib/sq in. 


NACA TN NO. 1279 


Fig, 7 



Fi9ure 7, — Comparison between tKeoretical and f inite— difference solution 
stresses in parallel-sided disk of 20-inch diameter rotating at 10^000 
rpm with temperature gradient that varies as fourth power of radius from 
600^ F at center to 1200^ F at rim. E, 30x10^ pounds per square inch; 
a, 0x10*^ (inches per inch) per V; \i, 0«3. 




